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KINETIC APPROACH TO SOLITON TURBULENCE 


C. M. Tchen 

City College of the City University of New York, New York, N.Y. 10031 


ABSTRACT 

The forced Schrodincer equ a tion for strong fluctuations describes 
the micro-hydrodymical state of soliton turbulence. It is valid for large- 
scaie turbulence where the internal waves can interact with the velocity 
fluctuations. It was derived by Tchen, as a fluid analogue, from the 
Navier-stokes equations for compressible flows. It is transformed into 
a master equation, to be subsequently decomposed into a macro-group, 
a micro-group, and a submicro-group, representative of the three trans- 
port processes of spectral evolution, transport property, and relaxation. 

The loss of memory in formulating the relaxation function gives the 
closure so that the transport property can attain its equilibrium. The 
kinetic equation for the macro-distribution is returned to the continuum 
by the method of moments to derive the equation of spectral evolution. 

The spectral flow is governed by the three transport functions as follows. HA 
transfer function from the (cubic) modulational nonlinearity forms 
a direct cascade. 2) A transfer function from the driving force, 

which enters as a convective nonlinearity in the phase-space, produces 



a reverse cascade for accumulating the solitons toward the large-scale 
end of the spectrum. 3) A coupling function represents the excitation 
of the solitons by the driving force which acts as a scattering (or 
emission) from velocity fluctuations. 

The following power laws 

E k 2 ~ k' m , N k ^ k" n 

are found for the soliton field intensity and density, respectively, 
with 

m = 0 , n = 0 

in the inertia subrange by reverse cascade, and 

m = 2 . n = 4 

in the coupling subrange at larger wavenumbers. These analytical 
predictions are verified by experiments. The forced Schrodinger equation 
is valid for fluids and plasmas. 
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1. Introduction: The Forced Schrodinger Equation for the Description 
of Soliton Turbulence 

The multi-scale motions in fluids and plasmas are more adequately 
described by the nonlinear Schrodinger equation* 

O t ' iv nV 2 + = 0 . (1) 

instead of the usual Navier-Stokes equation. The modulation relation 

<V 'Vi | V I 2. 

N = - A^Ie f 

makes the Schrodinger equation cubicaliy nonlinear. Here E is the 

/Vi 

longitudinal field and is related to the density fluctuation N by (2). 
The density is normalized such that 

N = 1 . (3) 

The over- bar denotes an ensemble average, () denotes a fluctuation, 
and A is the operator, the constants are 

_2 

having the dimensions of viscosity, frequency, and (acceleration) 
respectively. 
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The representation of the many-scale motions by the Schrodinger 
equation is well recognized in plasmas, because of the presence of the 
ion and electron species, and therefore the Schrodinger equation was 
first derived for plasmas. The analogous representation for fluids, called 
the fluid analogue, is not evident. In fluids of large-scale motions 
where the compressibility is not negligible, the longitudinal density 
waves may couple with the transverse velocity fluctuations. The internal 
gravity waves of large scales may also couple with the velocity fluctuations 
and transmit the effects to smaller-scale turbulence. These fluid motions 
of many scales justify the fluid analogue of the Schrodinger equation 
and its derivation from the Navier-Stokes equation for compressible 
fluids. The reason for having a fluid analogue lies in the fact that 
the Navier-Stokes equation is formally tied down to small scales by 
the nonlinearity u • V u of the gradient type, where u is the fluid 
velocity . The Schrbdinger equation has the cubic nonlinearity 

i -v * 2 ** 

V E I E not of the gradient type, and is found more suitable 
for the description of large scales. 

An examination of (1) reveals that the modulational nonlinearity 
yields a direct cascade, and should require a sink for balancing the 
spectral flow. An artificial sink would introduce an ambiguity. Hence 
the Schrodinger equation in the original form (1) cannot describe turbulence. 
The reason for this difficulty lies in that the original derivation was 
intended for weak fluctuations. The fluid analogue for strong turbulence 
requires a new derivation. 

We have derived a generalized Schrodinger equation in the form 

3 
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as follows: 


( a t - « y n 7 2 + tl<JA N)E = -IX , 


with a driving force X , that obeys the Poisson equation 


(4) 


V X = -icj r , (5) 

2 n 

and contains the scattering (or emission) 

r = A 77 uu , (6) 

by the velocity fluctuation u . The equation (4) will be called the 
"forced Schrodinger equation". We consider isotropy by assuming 

1=0, X = 0 . (7) 

In conclusion, the forced schrodinger equation, but not the standard 
Schrodinger equation and the Zakharov equations, will develop 


turbulence. 



2. Master Equation for Soiitons 

For the evolution of the spectral distribution of field fluctuations, 
the knowledge of the eddy transport properties is necessary, for which 
a kinetic theory is most appropriate, like in any transport theories 
of gases and plasmas. A continuum theory of correlations, as based 
upon the forced Schrodinger equation, would immediately be faced with 
the difficulty of high-order correlations. The most well known kinetic 
msthod for treating collective phenomena is the statistial method of many 
bodies by Bogoliubov. It introduces a master equation of N particles 
that interact with a field , and by successive integrations, 
generates a hierarchy of equations for the many-particle distribution 
functions, called the BBGKY hierarchy. For the spectraltevoiution, 
the closure that leads to the kinetic equation for the pair-distribution 
is necessary. This is considered to be a difficult task in statistical 
mechanics. In addition, the field of interaction E(t,x) that is a simple 
Coulomb field in plasmas becomes much more complicate in soiitons, 
as having to satisfy the nonlinear partial differential equation (4) in 
view of its self-consistency. We also wish to avoid the difficulty from 
the kinetic equation of pair-distribution function. 

For these reasons, we devise a group-kinetic method. The 
forced Schrodinger equation (4) can be considered as the first moment 
of the following master equation 

O t + A L)f(t,x,E) = 0 . (8) 

The differential operator 
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S * l n * 4 


(9) 


L = 


consists of three conponents 


L v ' L X=- iX '' 3 • 1=^1 


The function 


( 10 ) 


f(t,x,E) = T [e - E(t,x)] 


( 11 ) 


in the phase space has E as an independent variable. The zeroth 
moment gives 


/ dE f = t, jTdE f = 1, y"dE 7 = 0, 


(12a) 


and the first moment gives 


I 


dE 


A 

E f 


= E, E= fdE E f , E =[ dE E f 

«Vi / «■» V MM* 


(12b) 


The master equation in the form (8) describes the microdynamical 
state of turbulence. This can be shown by reverting it into the forced 
Schrodinger equation through the moment operation. 

Since the field anplitude E is complex, the distribution function 

** /V A* 

f(t,xE) will be complex too, while the functions N and X are real. 
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3. Kinetic Equation of Soliton Turbulence 


For soliton turbulence, we need to derive a kinetic equation 
that explicitly shows the interaction between the large and small scales. 
Since both the forced Schrodinger equation and the master equation 
(8) as a kinetic correspondent describe the microdynamical state of 
soliton turbulence in too many details that are unnecessary for a statistical 
treatment, a procedure of coarse-graining is needed. The decomposition 
into Fourier components will not suffice. In compliance with the separate 
statistical roles, we decompose the fluctuation 

0 = 0° + O' (13a) 

o 

into a macro-group (} and a micro-group ()', and by re-scaling 
the micro-group 


O' = O 1 + ()" (13b) 

into a submacro-group ()^ and a submicro-group ()" . The three 

groups with operators 

A° , A 1 , A" (14) 

represent the three transport processes of spectral evolution, transport 
property, and relaxation. The loss of memory in the relaxation makes 
the transport coefficients approach their equilibrium, and yields the 
closure. The three groups are separated by the three durations of 

1-8 
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correlations 



CO CM 


? > T 


(15) 


The superscripts 


o * 1 




(16) 


with square brackets denote the deterministic statistical properties 

as derived from the fluctuations (13) and from ensemble averages. 

° 1 

We scale the master equation by the operators A and A to 
obtain the system of equations 


o * ft 7 o o — 

( + A L -£ U if = - L f 

and 

(> t + A* L -Z )f = - L'(f + f ) . 


07) 


08) 


The eddy collisions are: 

If/ 11 f° - - A° lV 

f 1 = - A 1 L"f" . 


(19a) 


(19b) 


1 
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The submicro-distribution 


f" 




( 20 ) 


is a cluster of high-order groups that represents the frictional environment 

i f| 1 f«l 

in which f evolves. The collision coefficients and may 

act as integral operators. 

The composition of is determined by (18). Upon integrating, 
we have the distribution function 


/. 


dr /V(t,t-T JL^t-T) 


(f+f°) 


t-r 


(21a) 


from which we calculate the eddy collision 



lV 



t 

4* LV) A Uut-T) L.\t-T)(f+f 0 ) 

X-T’ 


(21b) 


Since the groups 


1 0 
A 1 and A 


have the time scales 

rn c°i 

T c and \ 


by (15), we can replace the upper limit t of integration 



by oo without loss of generality, and obtain an asymptotic collision 
coefficient 



<V(t) A^t.t-T) L^t-T ))> 


(22a) 


that is deterministic. Here we have written 
lJ(t) = L*(t,x,E) . 

for the sake of simplicity, and A 1 is the operator of evolution. In 
an analogous manner, we define 

[A i-*oa 

5 = J dr^L'(t) Aht't-V) L'(t-T)) . (22b ) 

0 

By noting 

- A* LI' = - aVl 1 +L (2 >L (2) 


we can show the equivalence between (22a) and (22b), and write 


e 


cu 



(23) 


J 
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With the collision coefficient thus determined in the form (22b), 


the equation of evolution of the macro-distribution (17) becomes 


< l 


t 


+ 


O A o 

A L)f 


L f 



( 24 ) 


o 

It is explicit in f , and is called the kinetic equation of turbulence. 


4. Schrodinger Equation of Soliton Turbulence 


By taking the moment of the kinetic equation (24), we revert 
to the continuum and derive the following Schrodinger equation of 
soliton turbulence: 


7 o o o o o 

( 2>t - 1 ^n V + 'z^n A N )E = ‘ iX + J ' 


(25) 


with 



(26a) 


or equivalently 


o 

J 



o 

f 


( 26 b) 


by (23). It has ail the constituents of the original forced Schrodinger 


j 
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equation (4), and, in addition, the eddy stress (26). The latter 
represents the statistical effects of the colliding micro-eddies. 

By writing 



(27) 


into two components as corresponding to L^ 1 and L^‘ , we transform 
the eddy stress (26b) into two parts as follows: 


o o o 

J = + "*X ' (28a) 


with 

o r ^ f»l o o t so 0) o 

J N = JdE E *? N f , J x = J dE E € x f . (28b) 

The limits of integration are understood to extend from -o» to 00 . 

The coefficients of collision are integral operators, and the operator 
symbol is understood. 


5. Spectral Flow 

o* 

Upon multiplying (25) by E and adding the complex conjugate 
part, we obtain the energy equation for the solitons, as follows: 


j 
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X 

L 




T 


to! 





( 29 ) 


Homogeneity and isotropy are assumed. 

The transport functions are written in the order of their importance 
with increasing wavenumbers. 

The coupling function 

w = -L [ - i/X • E > + i < X E > j (30) 

excites the solitons by the driving force. The transfer function 

foj ft O o* I V 

T x = -(l J x • E V {31a) 

for the reverse cascade accumulates the solitons toward the large-scale 
end of the spectrum, and the transfer function 

T N f ° ? O) (31b) 

for the direct cascade causes a disintegration of the large eddies 
into smaller eddies. The vertical bars represent the absolute values. 

The remaining convective terms 


3 
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(32) 


( - i v> n V 1 + i £v n A° N°)E° 


in the left hand side of (25) do not contribute to the spectrai flow, 
on account of their imaginary character. 


6. Theory of Eddy Transport 


6.1. Eddy aiffusivities 


It will be convenient to introduce the following eddy diffusivity 
operators 



for transforming the collision coefficients into 



(33a) 

(33b) 

(34) 

(35a) 

(35b) 
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and the transfer functions into 


Vs-JdE E.^<|e‘Y)> 

= (i"n) 2 JdE E.K N r, ^|E°Yl) 

T x Co, = -/ dE . !' ? x C ' 1 <i? OY, ^ > 

= f dE E3K^-i<^S°*f°|) ■ 


(36a) 


(36b) 


Use is made of (28), (31) and £35). The approximation retains those 
components of diffusivity (33) originating from the auto-correlations, 
and neglects others from cross-correlations. 

The diffusivities upon which the transfer functions depend 
are the T -integrations of the Lagrangian correlations, the orbit being 
characterized by the fluctuating evolution operator /V . For the sake 
of simplicity and as part of our attempt of memory loss, we approximate 
the operator by its average 

A' *A 61 , 

and write the diffusivities (33) in the Fourier form, as follows 


i 
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(37a) 



Here 

>5 (T /M>" 


is the coefficient of Fourier truncation in three dimensions within a 
length interval 2M which may be as large as desired, and 


N 


X 


[•1 2 

k 

[0 2 

k 


J dk" > ^ N'(k") N'(-k")^> 

jfdk" ^ <^X , (k”)-X , (-k")) 


are the spectral intensities of 


C»1 
k=0 - 


N- and X- fluctuations, such that 



The intensities (38) serve as operators. 


j 


(37b) 


(38a) 

(38b) 
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The submicro-time 




M 

h CC) 


(39) 


is the relaxation time for the approach of the diffusivities (37) to 
equilibrium. The same relaxation holds for both diffusivities^ since 
they share the identical operator of evolution /\* • The orbit function 
is 

[«"] Oil C"1 

h L (T) = h E (T7) h N (T ) h x u (x ) 

0»1 

= h E (T) h N (x) (40) 

with the components: 


h (T) = exp(- J-ik"- E T 2 ) (41a) 

c 4 w 

C"l y* f«l 

h N (T ) = exp . (41b) 


The first component relates to a streaming by E in weak turbulence, 

r f**J 

and the second component relates to an internal modulation by ^ < 0 

in strong turbulence. The third component 



= exp 



T 


is omitted in view of its external and divergent character. 



The orbit function is the result of integrqtion of the equation 
(18) for f 1 by the use of the evolution operqtor /J that entails 



From the characteric equations of the partial differential master 
equation, one can determine the exact orbit with more orbital components 
than those listed in (41a) and (41b). We have selected the major ones 
that represent streaming and modulation only. With the assumption 
of similarity between 


£ 


N 




and 


s*> M 
G N 


(42) 


a loop is formed, causing a loss of memory and yielding a closure. 

Upon substituting (40) into (39), and subsequently into (37), 
(35) and (36), we transform the diffusivities into 



C" 1 ) 

h E ( t) h N ( t ) j- 

[« ] 'j 

h E i ~ C) h N (T) f ' 


and the transfer functions into 



(43a) 


(43b) 


(44a) 


(44b) 


4 
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where 


H 

®k < 


r > ,/« . 


f.J o* o \ 

•h £ (-c)h N ( T ) E°*f°r> (44c) 


is the modified energy, as obtained by weighting the field energy 




by the orbit function. 


By noting that 


2. % _ k» 2 - 4 


^ h £ (T ) = - 


hpfT ) 


(45a) 


and by partial integration, we can write 


f " 1 \ Ic"' 1 ' n L" J 

3-h E (T)h N L J (x)c) =-JL T 4 h £ (r) h N (r) (45b) 


,.2 


C--1 


and 


t-*CD 

' C«3 

dr ^ -h E (r ) h N c> 


..,2 r t -»« 


r 4 [«] 

I dx x h E (-c )h^ ( x ) 

■Jn 


(45c) 


The approximation of keeping the highest moment in T has been used 
in (44a) and (44b). 


4 
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The change of signs in the two transfer functions (44a) and 

[ 6 ] r°] 

(44b) indicates that is a direct cascade and that is 

a reverse cascade. 


6.2. Fluidization as a means of closure 

We recall that the decomposition into the three groups with 
operators 

A° , A* , A" 


corresponds to a decomposition into three transport processes of evolution, 

diffusivity, and relaxation. Firstly, the macro-scaling yields the spectral 

[07 

evolution (29) with its coupling function W and the transfer functions 


Co] 


C°1 


N 


. By definition (30), the coupling function governs 


the spectral flow between the two macro-groups X and E .By 

[o] 

definitions (31a), the transfer function governs the spectral 

O 

flow from X 1 into E through the intermediary of the collision 

w 

>=> C'l 

coefficient b v . Similarly by definition (31b), the transfer function 

w 


T. 


N 


governs the spectral flow from E into N 1 through the 


Y> [i] 

intermediary of the collision coefficient b . Secondly, the collision 

C'i ro 

coefficients are proportional to the diffusivities and 

by (35). Thirdly, the approach to equilibrium and hence the closure 
depend on the time of relaxation , which enters into the diffusivities 

by (37). We have determined the relaxation by analyzing the orbit 

C'O 

function in two components h^(T ) h^ (x) fro m (40) and (41). 

They control the equilibrium of the diffusivities (43) and the trarsftr 
functions (44). 
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In strong turbulence the streaming by E is negligible, giving 

h E (T ) = 1 • (46a) 

The cluster of distribution functions (20) which shapes the frictional 
property of the medium in which f* evolves is now assumed to act 
like a fluid with a collision coefficient 


fill n 

^ [|,J 

C Nf ( H> = G N “'i -! 1 


(46b) 


E=0 


This means that the collision coefficient 


N 


E) in the kinetic 


representation with its individuality in E as an independent variable 


R 


will be reduced into a collision coefficient (t,x) in the fluid 


representation by losing the individuality in E . In this way, ^ 


C -7 


Nf 


ceases to be an operator. Hence the "fluidization” of the cluster 
of distributions yields the closure. 


6.3. Calculation of the collision coefficients, the diffusivities, and 
the transfer functions 


With the approximation of strong turbulence (46a) and the hypothesis 
of fluidization of the cluster (46b), we reduce (41b) and (44c) into 


N 


•n 


- exp £* 
Nf 


M 


(47) 


and 


i 
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(48) 


0. (T) = h c CC) h w (T) 


<1* 1 2 > 


At the same time we calculate the moments In ~0 


dr h P (r) h M CC) 


• ICI 


F d-c h E (r ) h N ^(T ) = c x 7 ~ 5 

0 t N X Nf 


with c Y = 4! , to transform (42), (35) and (44) into 


k n = N k i\t 1 > 


‘xV^ICT)- 


l 4 1 


Cl 1-51 


T t"l= -fc 1 ' 1 E^ 2 

’Nf Nf k 


_ l°l = - 'C C'l (°l 2 

T Xf Xf 

respectively. The subscript ( ) f denotes a fluidized property. 
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By writing the spectral intensity in terms of the integral of 


the spectral density, we have 


N, 


C ']2 




= 2 


dk" F N (k ,, J 


C'3 


X 


k 



F x (k") , 


( 53 ) 


where F N (k") and F x (k") are spectral densities of the N- and 
X- fluctuations. 

By the use of the notation (53), we can rewrite ( 51a ) in the 
form of the following integral equation 


2 / dk " F N (k,,) 


*C<J I * 1 

£ (k")j 

Nf I 


(54) 


. v=> Co 

for the determination of the collision coefficient (Z Kir . Here we have 

Nf 

Vo ll) 

assumed that U and (Z, have the same spectral structure 
Nf Nf 

except for the difference of arguments k and k" • The integral 
equation is solved to give 


St ' 1 

Nf " 


T"] 


u N k 
fZ k 


(55*) 


With this solution we calculate all the expressions in (50) and 
(51) to obtain: 


j 
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(55b) 



(56a) 


(56b) 


6.4. Coupling process 


The coupling function , as defined by (30), couples the 

o o 

two macro- fluctuations E and X .Its determination requires a 

**» *** 

fluid equation of macro-evolution in the form 



V?‘ 


i 

l n 


A N 


' ^Nf) 



(57) 


as obtained by (25) under the hypothesis of fluidization and by using 
the asymptotic value of the collision coefficient at large scales. 

Upon integrating (57), we get the macro- field 


E 



00 o 


O 


dT A (t,t- T ) X (t--e) , 


and calculate the coupling function 


(58 a) 


i 
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Here 


Co] 

K x is 


the trace of the diffusivity tensor, and 


(58b) 



( 59 ) 


is the evolution operator. An approximate estimate is given in the 
following. In the Fourier form, (58b) is 


X. W 2 
k 



dx 


h N (r) 


In the fluid approximation, we have 


( 60 ) 



( 61 ) 


and obtain 


W 


M 



y C°1 2 

A k 


( 62 ) 


j 
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7- Spectral Structure 


7.1. Equation of spectral balance 


Co] 

We have obtained above the transfer function for the 

ye f«7 

reverse Cascade (52b) with the coefficient of collision by (51b), 

[o] 

the coupling function W by (62), and finally the transfer function 

[°] T 1 1 

for the direct cascade (52a), with the coefficient by 


(55a) 


The three transport functions govern the spectral balance, as 


follows : 


£ 


E = 





_ C C' 1 F C°] 2 

■ Xf t k 



£•7 2 

k 


r t*7 c 

+ rof E k 


C°] 2 


(63) 


The function 


£ 


E 


2 > 


represents the instability of solitons. 

The inertia subrange by reverse cascade is governed by 


£ 

E 




(64) 


i 
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The coupling subrange at larger wavenumbers is governed by 


the spectral balance in the form: 


or 


f°l 


Co] 


w - t n = 0 , 


'Sv.J' 1 xJ-° -I I E.f 


Nf 


Nf 


= 0 


(65a) 


(65b) 


The asymptotic value of the collision coefficient is used in the present 

r°i 

subrange of larger wavenumbers. The amount of energy W that 

Col 

is produced by emission flows into the direct cascade T^ 

[ol 

Formally, a constant T^ would predict a secondary inertia 
subrange by direct cascade, if a constant dissipation could be found. 

Even so, the spectral cutoff by can occur before the direct cascade can 
be established. 


7.2. Spectral laws 


We take the spectral law 


[ ° 1 2 


X 


k 




( 66 ) 


for the driving force, where is the frequency scale . This 

spectrum is valid for a driving force that serves as a production of 

3 

internal density waves in compressible turbulence. 

The spectral structure in the inertia subrange by reverse 
cascade is governed by the spectral balance (64), under the conditions 

1-28 

j 



(2) and (66). We find the spectral laws 



«E< V’V V 5 


„ , S -V 6 -5 i -1,1/2 

C N ( £ E ^ X Hi ^ J 


C E 

C N 


.Vt 
A 




= c £ 2 



(67b) 


representing a flat maximum in the spectral plot. 

The spectral structure in the coupling subrange is governed 
by the spectral balance (65b). By again using (2) and (66), we find 
the spectral laws: 


E k 2 = FnI 2 ^x 6 k ‘ 2 ,68al 

N k 2 ^I^Nfl'^X 6 ’ 2 ^ • (68bi 


Here 




C'l 

Nf 


k=0 is an ' nterna * parameter defined by (55a). 
The spectral laws of intensities have their group notations omitted. 


8. Conclusions 

Although the equivalence between the Schrodinger equation and 

2 4—6 

the Navier-Stokes equations has been established, ' their nonlinear- 
lity differs. The former is cubically nonlinear and the latter has a 
nonlinearity of the gradient type. The Schrodinger equation is thus 
more suitable for the description of large-scale motions than do the 
Navier-Stokes equations. 


1 
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The nonlineqr Schrodinger equation and the more general 

Zakharov equations were originally derived for weak fluctuations in 

plasmas. 1 Their derivations neglected the scattering (emission) by 

velocity fluctuations. By generalizing to include finite fluctuations, 

we found the forced Schrodinger equation/ 

By the reverse cascade, the field energy is accumulated toward 

the large-scale end of the spectrum. The inertia subrange is followed 

by the coupling subrange, where the energy that is built up by emission 

is cascaded down toward the smaller eddies in a direct cascade. 

Consequently, the spectral intensity of field fluctuations falls with 

the power law k 2 , and the spectral intensity of density fluctuations 

-4 

falls with the power law k 

The role of iV is to cause a spectral cutoff of these power 
laws. The prediction of the power laws for the density intensity has 
been measured and is verified by True in the plasma experiments. 7 
See Fig. 1. The plot uses the spectral density F N (k) that is related 
to the spectral intensity by 


N 


C* 3 2 

k 


.co 


dk" F(k") 


( 69 ) 


The absence of the driving force from strong fluctuations 
will miss the mechanisms of accumulation and coupling. The transfer 
by direct cascade is the only mechanism which survives. On its own 
it cannot initiate a spectral flow and a spectral balance. This explains 
why it is so difficult for the standard Schrodinger equation to yield 

O 

a spectrum. 
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FLUID ANALOGUE OF THE SOLITON FORMALISM 


C.M. TCHEN 


City College of the City University of New York, New York, N.Y. 10031 


ABSTRACT 

In view of the many applications to turbulent phenomena in fluids, 
plasmas, optics, astrophysics and nerve systems, a large effort has been given 
to investigate the solitons analytically and numerically. The results 
of analysis by means of the Z-S equations (the nonlinear Zakharov equations 
and the Schrodinger equation) have not been encouraging, and have 
concluded that the solitons could not become turbulent and establish a 
broadened spectrum. An examination of the mathematical foundation of 
these nonlinear equations reveals that they were indeed intended for 
weak fluctuations and not for turbulence. By a nonlinear analysis that 
carries all the nonlinear fluctuations including those from velocities, 
we develop a general soliton formalism that emphasizes the dynamics 
of the two-scale motions for a parallel development of turbulence in fluids 
and plasmas. From the Navier-Stokes equations for plasmas and compressible 
fluids of two scales (fast and slow waves), we derive two equations of 
propagation of density waves. The fast wave is related to the fast field 
by the property of the spontaneous creation of the field by rarefaction, 

and the slow density wave is related to the field intensity by the property 
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of the ponderomotive force. These two properties, which have been taken 
for granted in plasmas in the Z-S equations, are not evident in fluids. 

We demonstrate them by a Lagrangien formalism and a kinetic method 
as consistent v/ith the Navier-Stokes dynamics. The first of the two properties 
transforms the hyperbolic equation of propagation of the fast density 
wave into a parabolic equation of evolution for the field-envelope as driven 
by the sound emission from finite velocity fluctuations. The second property 
gives the modulation. The emission and the modulation produce two 
nonlinearities. Our soliton system takes the generalized form of the Z-S 
equations by including the sound emission. It can be degenerated into 
the original Z-S equations by the omission of the sound emission in the 
weak fluctuation approximation. In our general formalism the solitons 
can be unstable and become turbulent, as occurring in the atmosphere. 

There, the measurments have shown that the solitons play a particularly 
important role in large scales, because the Navier-Stokes equations which 
are nonlinear in the form a gradient are not suitable. 


2-2 



1. Introduction 


There has been a great interest recently in the modulational 
instability and turbulence of solitons, as described by the Schrodinger 
equation 1 



( 1 ) 


with 


V 


N = 



( 2 ) 


where E g is the envelope of the field fluctuations, N is a density 
fluctuation, A is an operator of fluctuation, and 

' Ci ^ n ' and A 

are constant quantities having the dimensions of viscosity, frequency and 
_2 

(acceleration) 

The parabolic equation with cubic nonlinearity (1) has applications 
to a broad range of problems in nonlinear phenomena as connected 
with nonlinear optics, plasma, nerve systems, internal gravity waves, 

2-3 

j 



water waves and fluid dynamics. It has been proposed already in 1944 
2 

by Leontovich for analysing the nonlinear propagation of light and 

radio waves. The soliton description has become one of the standard method 

for treating nonlinear waves in multi-scale systems, where the fast waves 

are modulated noniinearly by slow waves. Its importance to instability 

of plasmas and to various problems of solid state has been demonstated 

by the large amount of analytical and numerical works in the current 

3 4 

literature. Its relevance to gravity wave and water waves ' , as well 
as to electronic excitations in long helical molecules ^ has also been 
demonstrated. 

The soliton description has been shown to be valid in a wide 
variety of other applications in fluids .rotating flows, acoustic turbulence, 
the dynamics of liquid sheets, thermal convective instability, Rayleigh 
instability, and instability of Poiseuille flow. In the atmosphere, many 
of these instabilities will develop turbulence. 

The importance of the cubic nonlinearity to fluid dynamical 
turbulence has been pointed out already in 1944 by Landau and the 
cubic parabolic equation in the real form has been proposed by Stewartson 
and Stuart 7 . 

The equivalence between the Navier-Stokes equations and the 

soliton formalism has been shown by several authors for fluids and plasmas. 

8 — 1 0 

The equivalence in fluid was based upon the Madelung 

transformation 




( 3 ) 


j 
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The derivation was formal, involving certain restrictive conditions 

(irrotational .flow and special barotropic property), and did not distinguish 

between scales. On the other hand, the equivalence in plasmas 1 did 

distinguish scales, but took the properties of space charge and ponderomotive 

11 12 

force as axioms, without going into their mathematical foundation 

The Schrodinger equation thus derived suffers from its basis from weak 

fluctuations. 

In the following, we develop a nonlinear theory of solitons to 
be valid for both fluid and plasma. The nonlinearities come from two 
sources : The modulational nonlinearity is based upon the ponderomotive 
force from the large-scale (or low frequency) density fluctuations, and 
the nonlinear emission of sound is produced by the strong velocity fluctuations 
of small-scales (or high frequencies). The density fluctuations of two 
scales act differently cr. the field : The density fluctuation of small- 
scale creates a spontaneous field-divergence, known as the space charge 
in plasma, and the density fluctuation of large-scale responds to the 
ponderomotive force, that is a known property in plasma and enters as 
the Madelung hypothesis (3) in fluid. We develop two theories to explain 
these two phenomena. 


j 
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2. Hydrodynamics of multi-scale motions. 


2.1. Scaling into slow and fast motions 

A 

A fluctuating quantity ( ) can be decomposed into an ensemble 

average ( } by the operator A and a fluctuation by the operator 

1 - A , where 1 is the unit operator. In multi-scale problems, we distinguish 
between a slow fluctuation ) and a fast fluctuation {* ) by the scaling 
operators 

A and A (4a) 


These operators add to 


I = A 




A 



(4b) 


and correspondingly, a fluctuating quantity has the following components 


<)= o * h 


(4c) 


A fast fluctuation 



( 5 ) 


j 
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has an amplitude, or envelope, that varies slowly. It is denoted by 

/v /v 


( ) and is obtained by the scaling operator A . The representative 

a 3 

frequency is (J . An infinitesimal rate of damping l written for 
convergence will be omitted eventually. The complex conjugate part 
is denoted by c.c. 

The velocity u in compressible fluid consists of a transverse 
mode and a longitudinal mode u . 


2.2. Basic equations 

We consider the equations of continuity and momentum as follows 


\ Si S'* 

K 'n- + V. m ol —0 

* * J 


( 6 ) 



+ V. n u. U- 

j i 


2 _ ^ 
■-C V. ol 

i 


+ 


(7) 


The longitudinal mode is represented by the collective field Ej. The 
fluid velocity 


U 


^ t 
u. ■+■ 


0 -r) 



( 8 ) 


i 
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may include both the transverse mode u* and the longitudinal mode 

/sP 

u * by putting cr = 0, or may represent the transverse mode only by 
putting c~ =1. The speed of sound is c in isothermal fluid, so 
that the pressure gradient is c v n . The density is normalized to 
unity as 



(9) 


We consider a homogeneous medium with . 


A 


u = o , 



no) 


It will be conveniant to transform (6) and (7) into the wave form by 

cross differentiation for eliminating the common term V ^ <x. 

d t J, ' 

yielding the equation of propagation 



r V n E. 
<} 4 



(ii) 


where 



( 12 ) 


j 
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is the scattering (or emission) function. 

Since the two forms of the wave equation (11), either with 
<7- z 1 or with V =0, are equivalent, we have the condition 


-Y'nE 
</ / 


A 

A. 


i 


defining the collective field 


E.. 

J 


The product 


(13) 


A A 


'Vi. £ 

J 


(14) 


is called the modulation. 

The scattering function 


* _ __ /V A ^ 

A, = VY: n. u, ia- 

~ — — 

* t ^ X 

- n + a 


(15) 


has two modes : the transverse mode is 



= VY (n U 


(16a) 


3 
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and the longitudinal mode is 


K 


2 


r A A A 

'n u u 



( 16 b) 


as denoted by the superscripts 

± $ 

U and ( ) , ( ' 7) 

respectively. 

By the second derivative W , the scattering function is 
controlled by small scales, so that we can approximate the total scattering 
by 


A. 


■ A 

W S u u • 


( 18 ) 


The longitudinal component has a fast fluctuation 


2c(j 

/is 717 


u 


6 a . 


( 19 ) 


j 
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with the amplitude 




( 20 ) 


and a slow fluctuation 


A/ 



77 : 




( 21 ) 


The index ( ) M denotes the slow fluctuation as distinct from ( )_. 

We assume that the correlations between different modes are negligible. 
Analogous definitions can be written for the transverse mode. Such an 
approximation is usual in sound propagation in a turbulent medium. 

The quantities 


n. 


r 

u, E, r 


( 22 ) 


A 

with ( ) refer to both the rapidly fluctuating motion and the slowly 

fluctuating motion , as distinguished by the notations 







and 


(23a) 






(23b) 
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with the characteristic frequencies 


CO and O (24) 

for the wave motions of high and low frequency scales, respectively. 

These two-scaled motions may be called the two species, indicated by 

a subscript { )^ 

r 

When the two scales (23a) and (23b) are treated as the two 
species, it might be convenient for a certain purpose to write the fields 

^ - i r S (25a) 

— ** *V w> 


i (25b) 

A/ A/ S A/ 


in terms of a common field IE , to be distinguished by the parameters 

£ and C in the ratio 

V. N 


! * » I 


(26) 


J 
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The coefficient 



has a + sign for a neutral fluid, and e^ = 1 , e = -1 
species of opposite charges. The parameters depend on the 
considered, e.g.on masses and charges in plasmas. 
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( 27 ) 

for two 
species 


4 



3. Hyperbolic and parabolic equations for the propagation of density 


waves. 


3.1. Notations 


For the fast species, we define 


/ \ '"V" ^ 

n - 1 + ^ ~r ^ f 


£T = E f- £ / With 



(28a) 


and 


A ^ * 



(28b) 


'V 

A 


A 

*- 


£ 


SC * X 

A E - E 


(28c) 


and stale the modulation as follows : 



(29a) 



(/ 






+ /) 


n fc 


(29b) 
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A -n. E 


^ t= 

- Tl E 


< r ZE> 


^ ~ A 

-- 77 E 


aZE — ^ - fi -f-A J £ & 




'Tt 


'V . y ^V 

E - /4 ^E 


(29c) 


(29d) 


in view of 


A + A = 1, for a slowly fluctuating function 
A + A + A = 1, for a rapidly fluctuating function 


and 


A ( ) = 0, in a homogeneous medium. 

Here 1 is the unit operator. 


(30a) 

(30b) 

HOc) 


3.2. Hyperbolic equation for the propagation of high-frequency 
density waves 

By applying the scaling operator A, we transform (11) with 
U~ - 1 into : 



The modulation takes the approximate form 


v 

** A 

A 'n E 
*■» 



■p /4 ^ 


3 
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( 32 ) 



from (29b), upon replacing n by N for coupling between the two 
scales of density waves. 

Upon substituting (32) into (31) we derive the wave equation 
for the fast mode, as follows ; 




A* 


A 


T u e driving force is 


(33) 


I / V A/ 

A = /lA t - 

* » t * * ~ « 

= /} ^ - V ^ . 

i „ * 

?//>t £ i ) 

= A + A J . 



(34) 


There we have neglected the compling n E. 
have made use of the definition (13). 


of different scales and 


3.3. Parabolic equation for the envelope of density waves 


The envelope of the fast wave has been defined by (5) and 

A/ 

is obtained by the operator A . Upon applying this operator and neglecting 
the second time derivative 'i ^ n g of the slowly fluctuating envelope, 
we transform the hyperbolic equation (33) into the parabolic equation 


j 
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I 


( 35 ) 



r. a. 




vNe^ --L(Z*S ) 

^ I V *■ cu ) , 


ZcJ 

>v 


with 


V 

-Tv 



and 


. •‘‘AC 

6 — U) < n v 


/tV 



3.4. Hyperbolic equation for the propagation of low 


density waves 


In an analogous manner and by applying A to (11), 
the wave equation for the slow species of frequency scale 







/ 


with 


Q 

M 



(36) 


(37) 

frequency 


we obtain 
<1^, as follows 

(38) 


(39a) 
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(39b) 


/V 




+■ 



A V 



/ 


by (13). 

On the other hand, the species of fast scale 
a density n that varies slowly according to the wave 





/V + A — 


J 4 


* 7" ? ! f* 


where we have written 





7nf 

' 4 


and 



from (13), (25) and (26). 




0 p may have 
equation 


(40) 


(41a) 


(41b) 


2-18 



Upon dividing (40) by , we have 


- 1 


■? ( 


V- 


z ) ~ 

/ " n = 


<P + 

>1 


< 1 - 


A 

•€ 

/v 


V E 


N 


(42) 


with 




~9 


(43a) 



(43b) 


By combining the wave equations (38) and (42) for the two 
slowly varying densities N and n = N , respectively , we obtain 






(44) 


with 



(45) 


Here we have neglected the term 





by 'f > 0, 


from (26). 
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is canceled in plasma by 


The term (1 + e n ) V • 

, but it survives in fluids by e n = +1 , from (27). 


4 
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4. Theory of ponderomotive force for plasmas and fluids 


The dynamical equation 



-ft 


(46) 


with 



/ 


describes the evolution of the longitudinal velocity under the collective 
field. It can be integrated to get 




' E(t-Z') 


(47) 


in the Lagrangian description. The Lagrangian field 


sfi-T-J = £[t-r,' fft-tj] (11 , 

is defined as the field at the time t -T along the trajectory. 
We can calculate the product 




l *2 

tA. 







(49) 
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The Fourier formula 



( 50 ) 


gives the integrand 

/» , , , pr M !,( A) b ')i 

E(t-T'jE(t-T'j - f -e 

JJ 

- l/cj't' + AJ T - ) 

X 

(51) 


In a homogeneous medium, we can take the average with respect to t 
in the interval of time 2T. It is assumed that the scale of t is larger 
than any characteristic time scales. We find 




T-$a> 



ifcj'i- cd’jt 

Z, 


y 


(52) 


The factor 

y = T/T 

is called the factor of Fourier truncation. After the time average, (51) 
is reduced into 


i 
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T~9CO 



B (t~xj E(t-x v J — fcfu) ' 

-v L-\. */ 


-Lcd'fc. -T ") 

■e 


X" E(tf) e(-<4') 


Z 2. -t'Qfc-T") ~i(T'-hZ") 

r - ^ -e. , (53) 


if the high-frequency field is peaked at £J n as its characteristic frequency, 
by (24), (5), and by writing 


Jk Efa’J i(-cS) = ff 


(54) 


Upon substituting (53) into (49), we obtain 


./-i 

u a — 


E e 



(54a) 


with 



-tCz’+z") 

C -e. 


(54b) 


The infinitesimal damping £ guarantees the convergence. The upper 
limit t of integration can be put equal to co without altering the value 
of the integral. 




2-23 



The double integral is now evaluated by writing 



_[ 


to find 



T T 

'K. -Tv 





I 

, , 2 . , 1 - 

oj +t 

'K 



This yields the tensor 


a) f ( __ 

Ol OL - 


-2 a 

CO EE 

■n 


( 55 ) 


(56) 


(57) 


to calculate the scattering functions 


A 

A 


t 


VV: 


* * 


4~ 2 vv : B E 


IQ « zr } cS 2 £ vv* be 

M ' n 

-/ -2. ~ >Z 

= l \ * f 

> *• 7 1 A.I 

5 A* 

- * l s , 


(58) 


(59) 


) 
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by (18) and (43a), finding the ponderomotive potential 


c = 





A 



In (58b), we have written 


( 60 ) 


O'* 


E 

L 


/V 



I rv <v 

— EE 



and 


E-E 

bw M 


I 

z 


A/ 



z. 


(61a) 


(61b) 


in d dimensions. The formula (60 ) for the ponderomotive force is derived 
for its general validity to plasmas and fluids. 

From (44), (43) and (59), it is seen that the density wave N 
is produced by an emission Q„ by the transverse component 
of the velocity fluctuations, and a source Q M by the longitudinal 
component of the velocity fluctuations in the form of a ponderomotive 
force. 


5 
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5. Spontaneous creation of the collective field by the density waves 


(a kinetic description} 


As a starting basis, we use the master equation 


(\ * C ) ?(% - 0 


( 62} 


for the distribution function 

Fft % v) ~ 'K *[v- , 

Wfc 


(63) 


with the differential operator 


l = + ifc*}- ■ l = 


(64) 


The kinetic description (62) replaces the dynamical description (46). 
By considering the operator S = 1 - A that gives the deviation 
from the average, we transform (62) into 


(2> t + )(F - - -* F • 


(65) 


By integrating, we obtain the distribution 



MM ^ 


( 66 ) 



and by taking the zeroth moment we find the density 



The evolution operator A governs the Lagrangian description of the 
field Se (t - T ) along the phase trajectory, such that 


= if [t-i, 


( 68 ) 


The position 


A 

X 




(69) 


along the trajectory consists of a portion x - V X from the streaming 

^ ~ A 

by v in the Euierian description, and another portion from the path J 
as depending on the field in the Lagrangian description. 

By the use of the Fourier expansion, the Lagrangian field (68) 

gives 

1 / L [ L ^ ^ 

** * ~ f Elcj'k) (70) 



3 
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The orbit function 


Ux ) - L (t) 

1 / v w E k 


has the following two components, by {69) . The kinematic component 
of streaming 



(72a) 


gives 


C>i Cc) = uk T A (t ) 

~ V ~ v 


(72b) 


The dynamic component 



0- 


(73) 


is characteristic of the oscillation and the scale of the field for the species^ 
If the characteristic frequency Cj^ dominates the dynamics of the oscillating 
field, the approximation is valid. The hypothesis that the Fourier component 
has a singular peak at without broadening is valid in our present 

derivation of soliton equation for the laminar flow, which subsequently 
may serve to describe the microdynamical state of turbulence. 
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Upon substituting (72b) and (73) into (71), and subsequently 
into (70) , we obtain . 


rj l(idt-k *) 

Afit-rJ jjd^cik < T l*k- aEp(id,k) 

= ^ z jj4ujk A uh • bEp(<J,k) 


-T 




v-Sfitx) 

^ / 


(74) 


The approximation of neglecting the streaming as compared to the dominant 
frequency has been made The characterization by Ct)^ specifies the field, 

r 

the distribution, and the density, as 

SEjj , ^ , ly 

From (67) , (74), we calculate the distribution function 


sr f ■- -j * ' 




(F+Fjv-h efo), (75) 


and the density 


r 




rt- 9 CO 

dr t 




( 76 ) 


j 


2-29 



by taking the moment. The integral calculated as 



- (Cu) +ih 
-e. “ 



( 77 ) 


reduces (76) into 






(78) 


By a change of variables from the specified field 


the common field 


SE 


by (25), we find 



into 


vie =- y 




(79) 


for one single species, or 


f-Se 

v** «%> 



(80) 


for many species. Here 


i 



X r 

* S/V 


( 81 ) 


i 
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is a factor independent of the species, by (26). 

For two species cf densities N and n and characteristic frequencies 
ci n and O n , (80) is 


V- 2 |E = ~ Y(X ^ + % ) • 

(82) 


In particular, we have 


V 


E 



for the fast field. 


(83a) 


and 


V 


v-E = - 



X t 





for the slow field. 


(83b) 


The relations (83) describes the spontaneous creation of the 
self-consistent field by the density fluctuations without going through 
a process of evolution. These relations are valid for fluids with 

A A a 

= 1 and for plasmas with e = -1 , e K , = +1 . 

A r n N 
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6. The forced Zakharov equations 


6 . 1 . Envelope equation for the field 

In (35) we have derived the envelope equation 
h . By the relationship (83a), rewritten as 

3 

t-, r 2 /x ' 

V E- - — co "h , 

4 . 'tu CL 

by the use of (25a) and (81), the envelope equation for 
formed into an envelope equation for E , with 

3 

v 

0 

from (37). The result of transformation is : 


f'V 


\7 y 2 '- -L 
** Z 


Cd f\ 


N) 


E 

Q. 


= A. 

2 'n, ci 


Note that, in view of 



in a homogeneous medium, we may omit the operator A 
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j 


for the density 


(84) 

n can be trans- 
a 


(85) 


(86) 


in (86). 



By integrating (86), we find 


( i 'a + vy 2 - foJA N )e = X 

v t 71 I '’v / *2 /v -i 


(87) 


The forcing function is X , such that 


v-x,= 


-fOJ A , 

Z <X. 


(88a) 


or , upon integration. 


xftxj= 


U y 1 


4-7 T J " Jx- X'j 


A/ ✓ 


(88b) 


In the absence of the driving force, i.e X =0, (87) is 

u 

degenerated into the first equation of the Zakharov equation 


fi) 4 W)e = 0 (89) 

v t "*■ 2. A ' 

With the condition (2), the equation in the form (89) is also called the 
Schrodinger equation, as written in (1), and the equation in the form 
(87) will be called the forced Schrodinger equation. 

6.2. The second equation cf the Zakharov system with the driving force 

The density N is governed by the wave equation (44). By 
the use of the Poisson equation (83b), rewritten as 
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i 




: ~ 0 + v<) . 

by (25) and (81), we transform (44) into 

+ ( , + - f W)(< + ^/-Cj ^ A/ 

s ^Q M = v*g 

- A ej V 2 A ( E a | (9i) 

«*s» 

by writing, from (60) . 

t = ^Si 4 r *. (92) 

Vta 

with 

a -N t &*£)"' 

The approximation in (91) is made on the basis of the ur.equaiities 

Q <*< A 

71 

and 

O, « CO 

bJ ^ 

that are small by the factor £ « 1. 


J 
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7. Conclusions 


The hydrodynamics of solitons is governed by the system of 

a/ 

equations (87) and (91) for the variables E g and N , respectively. 

The system contains the modulational nonlinearity by the ponderomotive 
force at the large scale and the nonlinear emission of sound by the strong 
velocity fluctuations at the small scale. The derivation is based upon 
finite fluctuations. 

By the hypothesis of local nonlinearity that approximates (91) 
by (2), we reduce the system into the forced Schrodinger equation 

4 * I f J = X w 

^ H- 


The two nonlinearities can be clearly shown in the kinetic represt n- 
tation by the master equation 




(95) 


where the distribution function 




(96) 


has E as an independent random variable, and the differential operator 


f ~ AN - iX ^ 

L -v. 2 . _ 


(97) 


\ 
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contains two nonlinearities from the modulation by N and the sound 

<V 

emission by , as two convections in the phase space. 

If we confine ourselves to weak velocity fluctuations, we will 

a/ 

miss the emission of sound by X g , degenerating the system of equations 
(87) and (91) into the Zakharov equations and the equation (94) into 
the standard Schrodinger equation 






vS8) 


Our system of equations (87) and (91) and our forced Schrodinger 
equation (94) are valid for both plasma and fluid turbulence, since the derivation 
uses a general argument. 
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